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When  the  first  satellite  images  of  the  earth  appeared 
in  the  early  1960's,  it  was  apparent  that  cloud  patterns 
associated  with  extratropical  wave  cyclones  have  remarkable 
similarities  during  comparable  stages  of  development. 

Proper  interpretation  of  these  cloud  patterns  allow  the 
location  of  frontal  zones  and  vorticity  centers,  as  well  as, 
the  stage  of  storm  development  to  be  determined. 

Furthermore  these  patterns  agreed  extremely  well  with  the 
classical  models  of  clouds  and  weather  associated  with  wave 
cyclone  life  cycles. 

Although  synoptic  scale  cloud  patterns  over  the  entire 
globe  have  been  monitored  daily  by  meteorological 
satellites,  their  use  in  numerical  weather  analysis  has  been 
limited  by  the  difficult  task  of  making  quantitative 
inference  about  the  mass  structure  of  the  atmosphere  from 
cloud  field  data  (Broderick,  1969)*  Prerequisite  to  the  use 
of  satellite  observed  cloud  field  data  in  numerical  weather 
analysis  is  the  establishment  of  consistent  relationships 
between  cloud  patterns  and  parameters  which  are  readily 
transformable  into  data  suitable  for  a  numerical  analysis 


scheme . 


1 . 2  Objective 

The  primary  purpose  of  the  present  study  is  to 
formulate  an  objective  statistical  technique  to  derive  300 
mb  geopotential  heights  and  1000/300  mb  thicknesses  from  a 
combination  of  satellite  observed  extratropical  spiral  cloud 
patterns  and  parameters  computed  from  the  initial  height  and 
thickness  fields. 

1 . 3  Approach 

An  investigation  was  undertaken  to  establish  the 
statistical  relationship  between  the  negative  geopotential 
extrema  associated  with  short  (synoptic  scale)  wave 
disturbances  and  parameters  derived  from  satellite  observed 
spiral  cloud  patterns  and  the  planetary  scale  long  wave 
(smoothed)  geopotential  field. 

The  approach  to  the  problem  consisted  of  the  following 
steps : 

i)  a  detailed  review  of  previous  studies  to  aid  in 
developing  methodology, 

ii)  development  of  methodology  and  satellite  analysis 
scheme  based  upon  past  studies  and  present  technology  and, 

iii)  development  and  verification  of  the  multiple 
linear  regression  equations  for  estimating  the  magnitude  of 
the  negative  extrema  in  the  height  and  thickness  fields. 
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Chapter  2 


Background 

2 . 1  Satellite  Data  for  Numerical  Weather  Prediction 

During  the  first  decade  of  satellite  meteorology  (1960- 
70),  the  operational  value  of  satellite  imagery  in  synoptic 
analysis  was  demonstrated.  In  the  same  period,  weather 
forecasting  was  moving  into  the  age  of  numerical  prediction 
where  the  impact  of  satellite  imagery  was  limited  due  to  the 
difficulty  of  converting  imagery  into  information  useful  to 
numerical  analysis  (Bizzari,  1982).  Thus,  the  use  of 
satellite  imagery  has  been  limited  primarily  to  short  range 
forecasting . 

The  second  and  third  decades  of  satellite  meteorology 
have  seen  considerable  efforts  in  the  development  of 
technologies  and  techniques  for  providing  quantitative 
satellite-derived  data  for  numerical  weather  prediction 
models.  For  example,  the  input  of  satellite-derived 
vertical  temperature  profiles  and  cloud  tracked  winds  in 
numerical  prediction  models  have  had  a  positive  impact  on 
numerical  forecasts  (Halem,  et  al,  1977;  Seaman  and  Hayden, 
1979;  and  Tracton,  1977).  One  of  the  problems  of  this 
quantification  is  that  the  satellite  data  tend  to  reduce  the 
amplitudes  of  the  synoptic  systems  (i.e.  versus  those 
depicted  by  radiosonde  data  alone,  Miller  and  Hayden,  1978;. 


The  input  of  satellite  data  into  numerical  analysis  of 
the  mass  field  based  upon  qualitative  information  from  the 
cloud  field  has  been  limited  to  date  by  the  dirficulty  in 
establishing  consistent  relationships  between  the  mass  and 
cloud  fields.  Previous  studies  relevant  to  this  problem 
have  dealt  predominantly  with  the  500  mb  level.  This  was 
the  primary  level  for  upper  level  analyses  during  the  early 
years  of  satellite  meteorology  and  provided  the  most 
complete  upper  level  analysis  to  establish  consistent 
relationships  with  the  satellite  observed  cloud  field. 

Early  approaches  to  this  problem  (McClain,  Broderick 
and  Ruzecki,  1965;  and  Bradley,  Hayden  and  Wiin-Nielsen , 
1966)  made  use  of  the  approximate  relationship  between  the 
vorticity  advection  term  in  the  quasi-geostrophic  "omega 
equation"  and  vertical  motions  at  the  500  mb  level 
associated  with  the  satellite  observed  cloud  features.  When 
a  given  synoptic  situation  was  reanalyzed,  the  Laplacian  of 
the  stream  function  field  was  then  altered  in  the  region 
influenced  by  the  satellite  observed  cloud  system  and  then 
melded  into  the  initial  Laplacian  field  neighboring  the 
region.  The  result  of  modifying  the  Laplacians  of  the 
stream  function  field  at  500  mb  based  upon  satellite  imagery 
is  not  known  until  the  actual  stream  function  field  is 
retrieved  by  solving  a  Poisson  equation.  A  problem  with 
this  method  is  that  the  final  stream  function  field  can  show 


modifications  outside  the  region  of  influence  defined  by  the 
cloud  features.  These  modifications  are  a  computational 
result  of  the  scheme  used  to  numerically  integrate  the 
Poisson  equation. 

Nagle  and  Hayden  (1971)  developed  a  quasi-objective 
technique  to  modify  the  500  mb  short  wavelength  component 
field  as  opposed  to  modifying  the  Laplacian  of  the  stream 
function  field.  The  short  wavelength  component  field  is 
derived  via  a  scale  separation  process  in  which  t  h° 
geopotential  field  is  separated  into  additive  short  and  long 
wavelength  component  fields,  i.e.  a  spatial  mean  flow  and 
superimposed  disturbances  (Holl,  1963)*  The  scale 
separation  method  does  not  require  the  numerical  integration 
of  the  Poisson  equation  and  avoids  the  problems  associated 
with  the  region  of  influence  because  of  the  additive 
properties  of  the  long  and  short  wavelength  fields. 

Nagle  and  Hayden  developed  regression  equations  from 
which  the  magnitude  of  the  negative  extrema  in  the  500  mb 
short  wavelength  component  field  can  be  computed  as  a 
function  of  parameters  derived  from  satellite  observed 
spiral  cloud  geometry  and  the  long  wavelength  component  cf 
the  500  mb  geopotential  field.  A  critical  assumption  in  the 
technique  is  that  the  long  wavelength  componen"  of  the 
500  mb  geopotential  field  can  be  accurately  resolved  by  “he 
current  observational  network. 


Results  indicate  that  the  magnitude  of  th- 


extrema  in  the  short  wavelength  component  field  can  be 
specified  with  a  proper  interpretation  of  the  spiral  cloud 
features.  Furthermore,  the  resulting  modified  500  mb  height 
field  can  improve  the  numerical  analysis  in  data  sparse 
regions . 


2 . 5  Extension  to  300  mb  and  1000/500  mb  Fields 

Recently,  increased  emphasis  is  being  placed  upon  the 
300  mb  level  as  the  primary  level  for  upper  air  analysis. 
This  is  a  result  of  the  numerous  commercial  and  military 
aircraft  which  routinely  report  wind  data  at  flight  level 
(near  300  mb)  with  the  aid  of  onboard  inertial  guidance 
systems.  These  data  often  provide  a  large  proportion  of  the 
real  data  over  the  oceanic  basins  for  use  in  the  numerical 
analysis . 

Well  developed  extratropical  cyclones  viewed  by 
satellites  are  depicted  as  broad  scale  spiral  shaped  cloud 
patterns  in  which  there  are  cirriform  clouds  that  are 
closely  linked  to  the  flow  pattern  near  the  jet  stream,  i.e. 
near  300  mb.  Therefore,  the  extension  of  the  previous  work 
to  the  300  mb  level  is  physically  sound. 

One  of  the  problems  with  inserting  bogus  height  data  at 
a  single  level  of  a  multi-level  numerical  forecast  model  is 
that  of  maintaining  vertical  consistency  among  the  various 
levels.  It  is  well  known  that  the  typically  deep,  large 
scale  vertical  motions  associated  with  major  extratropical 
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cyclones  can  be  partitioned  between  two  major  forcing 
functions:  The  vertical  variation  of  vorticity  advection 

and  the  Laplacian  of  thickness  advection  (Barr,  Lawrence  and 
Sanders,  1966;  and  Krishnamurti ,  1966).  Therefore,  the 

satellite  viewed  extratropical  cyclone  is  reflective  of  the 
flow  pattern  within  a  deep  layer  of  the  atmosphere.  It  is 
for  these  reasons  that  the  extension  of  the  previous  work  to 
the  1000/300  mb  layer  was  attempted. 

2 . 4  The  Scale  Separation  Technique 

The  scale  separation  technique  developed  by  Holl 
(1963),  as  an  objective  version  of  Fjortoft's  (1952)  manual 
technique,  provides  a  fundamental  representation  of  the 
subject  fields  by  a  superposition  of  specific  elemental 
fields  which  represent  indigenous  patterns.  Holl  defined  a 
new  representation  of  the  spectrum  of  a  field  and  developed 
a  technique  for  decomposing  the  field  into  additive 
component  fields. 

The  technique  successively  applies  a  smoothing  operator 
to  the  subject  field  until  the  amplitude  of  a  specified 
wavelength  component  is  reduced  to  some  percentage  of  its 
initial  value,  say  five  percent.  Appendix  A  contains  a 
detailed  discussion  of  the  smoothing  process.  Since  the 
smoothing  operator  is  linear,  all  subsequent  additions  and 
subtractions  yielding  components  of  the  original  field  are 


commutative . 


The  smoothing  filter  as  defined  by  Holl  is  given  by: 


1 LW 


=  z  +  c 
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V?z 


da  = 


Zo  "  Zsw 


(1  ) 


where  is  the  long  wavelength  field,  ZQ  is  the  original 

field,  Z  Sw  is  the  short  wavelength  field,  C  is  a  constant, 

and  a  is  the  smoothing  parameter. 

Holl  originally  developed  this  technique  for  use  with 

the  Northern  Hemisphere  Polar  Stereographic  grid 

configuration.  There  are  63  x  63  grid  points  for  a  total  of 

o 

3969  points.  The  grid  distance  is  381  km,  true  at  60  N.  Th 
map  factor  for  the  grid  is: 

9 

1  +■  sin  60 

m,.  =  -  ,  (2) 

J  1  +  sin  0 | j 

where  i> .  .  is  the  latitude  at  grid  point  i,  j.  This  grid 
configuration  was  used  for  the  present  study.  The  degree  of 
smoothing  used  for  this  study  was  a  =  5.0. 

The  amplitude  reduction  factor  as  a  function  of 
wavelength  (in  grid  lengths)  is  given  in  Figure  1.  The 
response  curve  is  that  for  first  order  smoothing.  Since  it 
is  desirable  to  have  a  much  steeper  response  curve,  third 
order  smoothing,  as  explained  in  Appendix  A,  was  used  for 
this  study. 

Examples  of  the  scale  separation  technique  applied  to  a 
300  mb  height  field  are  shown  in  Figures  2,  3  and  4.  Figure 
2  is  the  initial  height  field;  Figure  3  is  the  long 
wavelength  field  resulting  from  the  smoothing  process.  The 


Figure  1 . 


Amplitude  reduction  factor  as  a  function 
f  wavelength  (in  grid  lengths). 


The  short  wavelength  component  field  is  significant 
because  it  closely  resembles  the  relative  vorticity  field 
(Nagle,  et  al.,  1966).  The  spiral  cloud  systems  are  related 
to  the  short  wave  systems  moving  through  larger  scale 
patterns  (Whitney  and  Herman,  1968).  Thus,  satellite 
observed  spiral  cloud  pattern  geometry  can  be  used  to  modify 
the  short  wavelength  component  field.  The  adjusted  short 
wavelength  field  is  then  added  to  the  long  wavelength  field 
to  yield  the  modified  total  height  field. 

The  short  wavelength  component  defined  in  (1)  is  given 


by : 


'sw 


da . 
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From  the  calculus,  the  mean  value  theorem  for  integrals 

2 

states  that,  if  7  2  is  continuous  over  the  interval  [o,  a], 
then  there  exists  a  value  a*,  where  0  <  a*  <a  ,  such  that 
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where  Z(a*)  represents  the  partially  smoothed  field  at  the 


degree  of  smoothing,  a*. 


The  geostrophic  relative  vorticity,  ,  is  given  by 


cg  =  gf"1  ?2z  . 


(5) 


where  g  is  gravity  and  f  the  coriolis  parameter.  Therefore 
the  relationship  between  the  geostrophic  relative  vorticity 


and  the  short  wavelength  component  depends  upon  the 


2  2 

latitude  and  t Y\e  equality  of  V  Z  (a*)  and  7  Z0.  A  direct 


proportionality  between  the  magnitudes  of  these  fields  can¬ 


not  be  expected  because  of  the  wavelength  dependence  of  the 


Laplacian  operator,  although  the  patterns  of  the  V  Z(ct *) 
and  ^  fields  should  be  remarkably  similar.  The  location! 
of  the  field  extrema  and/or  zero  values  should  also  be 


similar  (see  Figures  4  and  5). 
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Chapter  3 
Data  and  Analysis 


3 . 1  300  mb  and  1000/300  mb  Field  Data 

The  field  data  were  provided  by  the  U.S.  Navy's  Fleet 
Numerical  Oceanography  Center  (FNOC).  The  height  and 
thickness  values  were  available  on  the  polar  stereographic 
(63  x  63)  grid  points.  The  scale  separation  technique  was 
applied,  with  a  =  5,  to  generate  the  short  wavelength  (SD) 
field  and  the  residual  (SR)  field.  The  planetary  vortex 
(SV)  field  is  obtained  by  applying  the  scale  separation 
technique  with  a  =  16.5  which  acts  as  an  extreme  low  pass 
filter.  The  long  wavelength  (SL)  field  is  then  determined 
using  the  following  relationship  with  the  original  field  Z, 


Z  =  SR  +  SD  =  SV  +  SL  +  SD. 


The  SR  field  is  usually  considered  the  long  wave  field  (see 
Figure  3)  but  is  actually  the  sum  of  the  SV  and  SL  fields 
computed  by  the  scale  separation  technique.  The  300  mb  SV 
field  corresponding  to  Figure  2,  3  and  4  is  shown  in  Figure 


The  following  parameters  were  computed  at  the  centers 
of  the  significant  SD  field  extrema:  The  latitude  and 
longitude,  the  height  value,  the  u  and  v  geostrophic 


velocity  components  and  the  Laplacian  using  a  five-point 
finite  difference  stencil.  These  parameters  were  computed 

for  the  SR,  SV  and  SL  component  fields  of  the  original  300 
mb  height  and  1000/300  mb  thickness  fields.  The  data  set 
consists  of  the  above  parameters  computed  for  both  the  1200 
GMT  analysis  and  the  12  hour  forecast  from  the  previous  0000 
GMT  forecast.  The  data  was  collected  on  a  daily  basis  for 
the  time  period  29  August,  1981,  to  19  September,  1982. 

Plots  of  both  the  300  mb  height  and  1000/300  mb  thickness  3D 
fields  for  the  1200  GMT  analysis  were  also  collected. 

3 . 2  The  Satellite  Data 

The  satellite  data  used  f<?r  this  study  was  acquired  by 
the  NOAA-7,  which  was  the  second  operational  satellite  in 
the  TIROS-N  series.  Launched  on  23  June,  1981,  the  NOAA-7 
satellite  is  a  sun-synchronous  (polar-orbiting) 
meteorological  satellite  with  an  orbital  altitude  of  850  kjn 
and  a  travel  speed  of  about  6.6  km/sec.  The  NOAA-7  Advanced 
Very  High  Resolution  Radiometer  (AVHRR)  is  a  five  channel 
instrument:  two  channels  in  the  visible  and  near  infrared 

(IR),  0-58  -  0.68  um  and  0.725  -  1.10  um;  one  in  the  IR 
window  region  at  4  um,  3*55  -  3-93  um  and  two  in  the  thermal 
IR  window  at  11  um,  10.5  -  11-3  um  and  11.5  -  12.5  um .  The 
imagery  used  for  this  study  was  from  the  visible,  near  IR, 
and  thermal  IR.  The  images  used  were  Northern  Hemisphere 
Polar  Stereographic  mosaics  composed  of  successive  satellite 
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passes  over  the  North  Atlantic  ranging  in  time  from  1300  - 
2000  GMT.  The  satellite  imagery  was  available  for  the  same 
period  as  the  field  data  collection  period  with  the 
exception  of  the  entire  month  of  March  1982.  The  images  for 
March  were  not  available  and  thus. not  included  in  this 
study. 


3 . 3  Satellite  Analysis  Scheme 

The  well  documented  correspondence  between  specific 
satellite  observed  spiral  cloud  features  and  synoptic 
features  associated  with  upper-tropospheric  relative 
vorticity  fields  was  summarized  by  Nagle  and  Hayden  as 
follows : 

1)  The  boundaries  of  cellular  convective  cloud  areas  and 
the  trailing  edges  of  frontal  cloud  bands  correspond  to  the 
zero  line  in  the  relative  vorticity  field. 

2)  The  locations  of  extratropical  spiral  cloud  centers 
associated  with  cold-core  systems  correspond  to  positive 
vorticity  maxima. 

3)  The  leading  edges  of  frontal  cloud  bands  correspond  to 
axes  of  maximum  of  negative  relative  vorticity. 

4)  The  apex  of  the  anticyclonic  curvature  of  the  leading 
edges  of  frontal  cloud  bands  correspond  to  the  negative 
relative  vorticity  maxima. 

These  relationships  permit  a  rather  complete 
specification  of  the  pattern  of  the  short  wavelength 


component  field  provided  a  proper  interpretation  of  the 
spiral  cloud  pattern  is  made-  The  magnitude  of  the  negative 
extrema  in  the  short  wavelength  component  (SD)  field  remains 
to  be  specified. 

The  spiral  cloud  patterns  were  analyzed  as  follows.  A 
quarter  hemisphere  polar  stereographic  NOAA-7  mosaic 

o 

(preferably  IR)  of  the  North  Atlantic  covering  from  10  E  to 

Q 

80  W  longitude  was  selected  corresponding  to  the  date  of  the 
field  data.  The  300  mb  and  1000/300  mfc>  SD  field  plots  were 
compared  with  the  satellite  image  to  identify  the  spiral 
cloud  systems  with  significant  and  distinct  negative  extrema 
in  the  SD  field.  If  a  spiral  system  was  evident  without  the 
presence  of  a  distinct  SD  feature,  the  SD  field  was 
considered  in  error  and  the  analysis  was  discontinued  for 
that  date.  This  was  an  effective  gross  error  checking 
procedure  which  prevents  spurious  parameters  from  entering 
the  statistical  analysis. 

The  latitude  and  longitude  of  the  center  of  a  spiral 
cloud  pattern  was  recorded  as  SLAT  and  SLON  to  the  nearest 
whole  degree.  The  latitude  and  longitude  of  the  center  of 
the  SD  extrema  as  computed  by  the  scale  separation  technique 
program  was  recorded  as  FLAT  and  FLON.  These  would  later  be 
used  as  identifiers  in  searching  the  field  data  file  to 
extract  the  other  computed  parameters  for  that  SD  feature. 
Since  the  satellite  images  are  mosaics  composed  of  successive 
NOAA-7  orbits  ranging  from  1300  to  2000  GMT,  a  systematic 


spatial  and  temporal  error  is  evident.  The  difference 
between  SLAT,  SLON  and  FLAT,  FLON  are  summarized  in  Tables 
1  and  2  for  the  300  mb  and  1000/300  mb  SD  fields 
respectively.  These  errors  are  also  partially  attributed  to 
errors  in  the  original  analysis. 

The  location  of  the  upper-tropospher ic  trough  (i.e., 
where  the  high  clouds  abruptly  end  or  have  a  distinguishable 
gap  along  the  major  frontal  cloud  band)  is  then  determined. 

A  line  is  then  drawn  from  the  spiral  center  to  the  point 
where  the  upper-tropospher ic  trough  intersects  the  outer 
edge  of  the  frontal  band.  This  distance  is  defined  as  the 
amplitude  (AMP)  of  the  satellite  observed  cloud  system. 

A  line  is  now  drawn  perpendicular  to  the  line  defining 
the  amplitude  and  passing  through  the  spiral  center.  The 
intersections  of  this  line  with  the  westernmost  and 
easternmost  edges  of  the  spiral  cloud  band  define  the 
wavelength  (RL)  of  the  satellite  observed  cloud  system.  In 
actuality  this  is  a  half-wavelength  when  considering  the 
ridge-trough  system  as  the  full  wavelength. 

All  distances  are  defined  in  terms  of  gridlengths  (381 

o 

km  at  60  N).  The  schematic  representation  of  the  above 
analysis  procedure  for  an  idealized  spiral  cloud  system  is 
shown  in  Figure  7. 

The  final  subjective  assessment  about  the  spiral  cloud 
system  is  the  stage  of  development  parameter  (1ST).  A 
modified  version  of  McClain  and  Broderick's  (1967)  stage  of 
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Figure  7.  Idealized  spiral  cloud  system 
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development  classification  system  was  used.  Figure  8  shows 
the  various  stages  of  development  used  for  the  study.  The 
stage  of  development  parameter  along  witlj  SLAT,  SLON,  AMP 
and  RL  were  obtained  for  each  spiral  cloud  system  that  could 
be  associated  with  a  specific  minimum  in  the  SD  field.  The 
total  number  of  cases  analyzed  in  the  above  manner  is  421 
for  the  300  mb  height  field  cases  and  420  1000/300  mb 
thickness  field  cases . 

3 . 4  Statistical1  Methodology 

A  regression  equation  of  the  form 

n 

y  =  b,x,  *  b2x2  +  ...  -b„xn  -  £  bjXj  (7) 

j-1 

is  used  to  apply  the  previously  discussed  parameters  to  the 
problem  of  specifying  the  negative  SD  field  extrema.  The 
special  case  for  which  x^  is  identically  one,  is  specified, 
which  means  that  b^  now  plays  the  role  of  the  regression 
constant  or  y-intercept.  Symbols  are  defined  as  follows*, 
y  =  negative  SD  field  extrema, 
bj  =  regression  constant, 
bj  =  regression  coefficients,  and 
Xj  =  potential  predictors  to  be  screened 
(  j  =  2 , 3 , . . . , n)  . 

The  stepwise  multiple  linear  regression  procedure  used 
in  the  current  study  was  the  Biomedical  Computer  Programs  P- 
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gure  8. 


Schematic  representation  of  the 
development  parameter . 


series  program  P2R  developed  at  the  Health  Science  Computing 
Facility  of  the  University  of  California  at  Los  Angeles.  A 
detailed  discussion  of  the  multiple  regression  technique  and 
program  P2R  are  given  in  Appendix  C. 

The  data  set  compiled  during  the  satellite  analysis 
procedure  had  to  be  matched  with  the  300  mb  and  1000/300  mb 
SD  data  sets  containing  the  various  long  wave  parameters. 

The  latitude  and  longitude  of  the  SD  extrema  (FLAT,  FLON) 
was  used  to  search  the  SD  data  file  which  also  contained 
FLAT,  FLON.  To  avoid  minor  differences  due  to  FLAT,  FLON 
being  listed  as  accurate  only  to  1/10  of  a  degree,  a  one 
degree  by  one  degree  latitude  and  longitude  box  was  used  to 
provide  some  leeway  in  the  searching  routine. 

The  300  mb  and  1000/300  mb  data  sets  were  now  complete. 
The  300  mb  data  set  contained  413  cases  while  the  1000/300 
mb  data  set  contained  373  cases,  each  with  the  date  time 
group,  the  observed  SD  value,  satellite  parameters  and  the 
long  wave  parameters.  The  data  sets  were  now  randomly  split 
into  dependent  and  independent  sets.  The  300  mb  and 
1000/300  mb  dependent  data  sets  contained  290  and  263  cases, 
respectively  while  the  independent  data  sets  contained  123 
and  115  cases,  respectively. 

When  analyzing  the  spiral  cloud  systems  it  became 
apparent  that  the  relationship  between  the  wavelength  (RL) 
and  the  SD  value  behaved  in  a  non-linear  fashion.  Based 
upon  heuristic  reasoning,  the  following  variables  with  non- 


linear  characteristics  were  added:  RL  *  SIN(SLAT), 

SQRT(RL),  RL/COS(SLAT)  and  RL  *  SIN  (SLAT).  Ratios  of  the 
amplitude  and  wavelength  were  also  added  as  suggested  by 
Lowe  (1983).  A  complete  list  of  the  variables  used  is  given 
in  Table  3. 

The  program  P2R  was  then  used  to  process  the  dependent 
data  sets.  The  variables  were  entered  or  removed  from  the 
regression  model  based  upon  the  F  method  (see  Appendix  B). 
The  F-to-enter  and  F-to-remove  limits  used  for  this  study 
are  4.0  and  3.9,  respectively  with  a  tolerance  of  0.05. 

Regression  runs  were  also  performed  for  300  mb  and 
1000/300  mb  data  sets  without  the  inclusion  of  the  satellite 
observed  spiral  cloud  parameters.  This  was  an  effective 
check  to  see  if  there  was  adequate  information  contained  in 
the  computed  long  wave  parameters  to  accurately  specify  the 
magnitude  of  the  negative  SD  extrema. 
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3.  POTENTIAL  PREDICTORS 

FLAT  LATITOOE  Of  THE  SO  EXTREMA 
f  LON  LONG  I TUOE  OF  THE  SO  EXTREMA 

SLAT  LATITUOE  Of  THE  CENTER  Of  THE  SATELLITE  OBSERVED  CLOUD  SYSTEM 
SLOW  LONG  I TUOE  Of  THE  CENTER  Of  THE  SATELLITE  OBSERVED  CLOUD  SYSTEM 

AMR  AMR L ITUOE  Of  THE  SATELLITE  OBSERVED  CLOUD  SYSTEM 

RL  WAVELENGTH  Of  THE  SATELLITE  OBSERVED  CLOUD  SYSTEM 

1ST  STAGE  Of  OEVELORMENT  i  1  -  FORMATIVE.  2  -  MATURE.  3  -  OCCLUDED 
SRU  U-COMRONENT  Of  THE  GEOSTRORHIC  WIND  DERIVED  FORM  THE  SR  FIELD 

AT  FLAT , FLON 

SR V  V-COMRONENT  Of  THE  GEOSTRORHIC  WIND  DERIVED  FROM  THE  SR  FIELD 

AT  FLAT. FLON 

SLH  GEOROTENTIAL  HEIGHT  Of  THE  SL  FIELO  AT  FLAT. FLON 

SLU  U-COMRONENT  OF  THE  GEOSTRORHIC  WIND  OERIVED  FROM  THE  SL  FIELD 

AT  FLAT. FLON 

SL V  V-COMRONENT  OF  THE  GEOSTRORHIC;  WIND  DERIVED  FROM  THE  SL  FIELD 

AT  FLAT. FLON 

SLL  F  I VE-ROI NT  FINITE  DIFFERENCE  ARRROX I  MAT  ION  FOR  THE  LAP  LAC  IAN 

OF  THE  SL  FIELO  AT  FLAT. FLON 

SVH  GEOROTENTIAL  HEIGHT  OF  THE  SV  FIELC  AT  FLAT. FLON 

SVU  U-COMRONENT  OF  THE  GEOSTRORHIC  WIND  DERIVED  FROM  THE  SV  FIELD 

AT  FLAT. FLON 

SVV  V-COMRONENT  OF  THE  GEOSTRORHIC  WINO  DERIVED  FROM  THE  SV  FIELD 

AT  FLAT. FLON 

SVL  FIVE-ROINT  FINITE  DIFFERENCE  APPROXIMATION  FOR  th£  LAPLACIAN 

OF  THE  SV  FIELO  AT  FLAT. FLON 
RLS  RL'SINCSLAT; 

RL  5  SQRTCRLI 

RLMA  RL»AMR 

AMOR  AMR/RL 

RLOA  RL/AM* 

RlDC  RL/COSCSLAT) 

RlS2  RL«SIN(SLAT)*SIN(SLAT) 
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The  basic  statistical  properties  of  the  variables 
included  in  the  290  cases  comprising  the  300  mb  dependent 
data  set  are  shown  in  Table  4. 

The  following  predictors  were  selected  for  use  in  the 
regression  model:  The  square  root  of  the  system  wavelength 
(RL5);  the  Laplacian  of  the  SL  field  (SLL);  the  height  of 
the  SL  field  (SLH);  the  Laplacian  of  the  SV  field  (SVL)  and 
the  U-components  of  the  SV  and  SR  fields  (SVU  and  SRU),  The 
coefficients  and  statistics  for  each  independent  variable  in 
the  equation  as  well  as  the  multiple  correlation 
coefficient,  the  squared  multiple  correlation  coefficient, 
and  the  adjusted  squared  multiple  correlation  coefficient 
are  shown  in  Table  5.  The  root-mean-squared  (RMS)  error  for 
the  observed  versus  predicted  SD  values  was  30.94  meters.  A 
plot  of  the  residuals  (observed-predicted)  for  this 
regression  model  is  given  in  Figure  9.  Upon  examining  the 
residual  plot,  it  appears  that  the  residuals  fall  in  a 
horizontal  band  centered  about  the  abscissa  without  any 
apparent  trends,  which  suggests  that  the  model  (nay  be  judged 
as  adequate  (Afifi  and  Azen,  1979).  The  probable  error  of 
the  correlation  coefficient  given  by 


VAR [ABLE 


MEAN 


STANDARD 

DEVIATION 


SKEWNESS 


KUR  TOS  I  S 


FLAT 

.82840 

. 15277 

- . 44665 

.62471 

FLON 

.  704  1  7 

. 35623 

.  1  7232 

-  1  .02394 

SLAT 

.84159 

. 15953 

-  .  47547 

.31218 

SLON 

.  50622 

.33495 

.18481 

-  .  92232 

AMP 

3. 39130 

1  .  27269 

1.01112 

1.547  1  2 

RL 

2.83716 

.  88580 

. 65806 

.  04054 

1ST 

2.35172 

. 48468 

.  52996 

-  1  .47690 

SD 

-217.96550 

53.35474 

- . 13033 

-  .  04560 

SR  U 

18. 74138 

8.36404 

.2781  1 

.42935 

SR  V 

-  .08276 

6 . 18397 

-.37603 

.32827 

SVL 

3.02100 

4 . 36309 

-  .87640 

1.21658 

SLH 

-80.23103 

91 . 74687 

-.50833 

.44344 

SLU 

5.76207 

7.31710 

-.38065 

1 . 65773 

SLV 

-  .05517 

5.71722 

-.42538 

. 36289 

SLL 

27.00623 

18.67461 

. 25430 

. 38237 

SVH 

85.83102 

189 . 75940 

-.41856 

-  .  43799 

SVU 

1  3  .  18621 

3.28316 

.43648 

- . 60756 

SVV 

.23448 

. 70480 

.22580 

1 .94812 

RLS 

2.08902 

.72151 

. 55330 

- . 18229 

RLS 

1  .66437 

. 25885 

.31383 

- . 40974 

RLMA 

10.52327 

7.18180 

1 . 80518 

4.77913 

AMOR 

1 .20176 

.27743 

1 .04864 

1  .  '0518 

RLOA 

.87356 

. 1 9004 

.51530 

1 . 02752 

RLDC 

4  .  43903 

1  . 56943 

.84057 

.  95840 

RLS2 

1 .5701  1 

.  63935 

.  43046 

- . 18950 

290  TOTAL  OBSERVATIONS 


TABLE  5 


5.  300  nb  dependent  data  set  regression  results. 


Variables  in  Equation 


VARIABLE 

COEFFICIENT 

STD.  ERROR 
OF  COEFF. 

STD.  REG. 
OOEFF. 

F-TO- 

REMOVE 

RL5 

-121 .63345 

8.2886 

-0 . 590 

215.37 

SLL 

-2.84518 

0.3077 

-0.996 

85.50 

SLH 

-0.58965 

0.075Q 

-1.014 

61.82 

SVL 

-7.42866 

0.9643 

-0.607 

59.35 

SVU 

4.73323 

0-7821 

0.291 

36.62 

SRU 

-0.82807 

0.3259 

-0.130 

6.46 

INTERCEPT 

=  -10.42469 

MULTIPLE  R  =  0.8146 

MULTIPLE  R-SQUARED  =  0.6637 


ADJUSTED  R-SCUARED 


0.6565 


where  R  is  the  correlation  coefficient  and  N  the  number  of 
cases  (Landsburg,  1964),  is  0.013. 

The  following  predictors  were  selected  for  the  no¬ 
satellite-parameters  regression  model:  The  Laplacian  of  the 
SL  field  (SLL);  the  height  of  the  SL  field  (SLH);  the  U- 
component  of  the  SL  field  (SLU);  the  latitude  of  the  SD 
extrema  (FLAT)  and  the  Laplacian  of  the  SV  field  (SVL) • 

Table  6  summarizes  the  regression  results  in  the  same  format 
as  Table  5.  The  RMS  error  for  the  observed  versus  predicted 
was  40.93  meters.  A  plot  of  the  residuals  for  this 
regression  model  is  given  in  Figure  10.  There  are  no 
apparent  trends  in  the  residuals,  but  the  scatter  about  the 
zero  line  is  large  enough  to  consider  the  model  inadequate. 
The  probable  error  of  the  regression  coefficient  is  0.023. 


The  basic  statistical  properties  of  the  variables 
included  in  the  263  cases  comprising  the  1000/300  mb 
dependent  data  set  are  given  in  Table  7. 

The  following  predictors  were  selected  for  use  in  the 
regression  model:  The  product  of  the  system  wavelength  and 
sine  of  the  system  latitude  (RLS);  the  Laplacian  of  the  SL 
field  (SLL);  the  stage  of  development  parameter  (1ST);  the 


TABLE  6 


300  mb  dependent  data  set  regression  results  (no 
parameters ) 


Variables  in  Equation 


STD.  ERROR 

STD.  REG. 

F-TO- 

VARIABLE 

COEFFICIENT 

OF  COEFF. 

COEFF. 

REMOVE 

SLL 

-4.32998 

0.3684 

-1 .515 

128.14 

SLH 

-0.82326 

0.0940 

-1.416 

76.64 

SLU 

-1.81360 

0.4362 

-0.249 

17.28 

FLAT 

-142.29917 

35.5345 

-0.407 

16.04 

SVL 

-5.10234 

1 . 6776 

-0.417 

9 . 25 

INTERCEPT 

=  -23.33036 

MULTIPE  R  =  0.6414 

MULTIPLE  R-SQUARED  =  0.4115 


ADJUSTED  R-SQUARED 


0.4011 


RESIQUflL  { OBSERVED-PREDICTED  3 


VAR [ABLE 


MEAN 


STANDARD 

DEVIATION 


SK-E-NESS 


kur res  I  s 


FLAT 

. 91340 

.14645 

- . 13985 

.54160 

FLON 

.73611 

.  39217 

.08333 

-  [  .042  19 

SLAT 

.83976 

.15368 

-  .  22732 

.  32350 

SLON 

.61395 

.35139 

. 1 9790 

- . 74965 

AMP 

3.55163 

1  .  32945 

. 86974 

.81578 

RL 

2 .93490 

.94653 

. 74317 

.1  4643 

1ST 

2.31939 

.47433 

.  66749 

-  1  . 26407 

SO 

-176.23950 

55 . 30478 

-  .  22924 

.  85936 

SR  U 

17.23193 

6.50968 

. 66252 

. 38902 

SR  V 

-.97719 

5.83091 

- . 91658 

1  .  10230 

SVL 

3  .  13065 

3 . 78742 

-  .59761 

.  67441 

SLH 

-50.23193 

78 . 36052 

-  .68650 

1 .87620 

Slu 

4.05323 

5.35694 

. 33100 

.  8908  1 

SL  V 

-  .  82890 

5.23563 

- . 80441 

1 .01086 

SLL 

19.76711 

1 7 . 20300 

. 26900 

1 .02785 

SVH 

64.71863 

184 . 98930 

- . 42864 

-  .  507 1 8 

SVU 

13.36882 

3  .  19887 

.06318 

-  .  49276 

Svv 

.  16350 

. ’7990 

- .  1  9728 

1  .  82237 

RLS 

2.  15556 

.75189 

. 7  1720 

.  16032 

R  L  5 

1.69166 

.27048 

.37551 

-  .  27238 

R  L  MA 

1  1  . 46225 

7.92125 

1  . 59687 

2  .  99606 

AMOR 

1  .21578 

. 26360 

.82531 

1 . 33291 

RLOA 

. 8604 1 

.  1  848  1 

.  76750 

1  . 6036  1 

RLOC 

4 .57093 

1  .  64756 

. 90890 

. 96800 

RLS2 

1.61331 

. 65670 

.61915 

.12310 

263  TOTAL  OBSERVATIONS 


height  of  the  SL  field  (SLH);  the  Laplacian  of  the  SV  field 
(SVL)  and  the  U-components  of  the  SV  and  SL  fields  (SVU  and 

SLU).  Table  8  summarizes  the  regression  results  in  the  same 
format  as  Table  5,  the  RMS  error  for  the  observed  versus 
predicted  was  36.02  meters.  A  plot  of  the  residuals  for 
this  regression  model  is  given  in  Figure  11.  There  are  no 
apparent  trends  in  the  residuals,  and  the  scatter  appears 
acceptable,  which  suggests  the  model  is  adequate.  The 
probable  error  of  the  correlation  coefficient  is  0.018. 

The  following  predictors  were  selected  for  the  no¬ 
satellite-parameters  regression  model:  The  heights  of  the 
SL  and  SV  fields  (SLH  and  SVH),  and  the  U-component  of  the 
SR  field.  Table  9  summarizes  the  regression  results  in  the 
same  format  as  Table  5.  The  RMS  error  for  the  observed 
versus  predicted  was  48.51  meters.  A  plot  of  the  residuals 
for  this  regression  model  is  given  in  Figure  12.  There  are 
no  apparent  trends  in  the  residuals,  but  the  scatter  about 
the  zero  line  is  large  enough  to  consider  the  model 
inadequate.  The  probable  error  of  the  regression 
coefficient  is  0.032. 


STD.  ERROR  STD.  REG.  F-TO- 
VARIABLE  COEFFICIENT  OF  COEFF .  COEFF.  REMOVE 


RLS 

-33 .62609 

3.8860 

-0.457 

74.88 

SLL 

-3.22387 

0.4249 

-1 .003 

57.57 

1ST 

38.99478 

5.6850 

0.334 

47 . 05 

SLH 

-0.51486 

0.1019 

-0.729 

25.52 

SVU 

4.07143 

0.8280 

0.235 

24.18 

SLU 

-2 . 16366 

0.5317 

-0.210 

16.65 

SVL 

-2 . 38556 

1.1721 

-0.163 

4. 14 

INTERCEPT 

=  -194.5239 

MULTIPLE  R  =  0.7588 

MULTIPLE  R-SQUARED  =  0.5758 


ADJUSTED  R-SQUARED 


0.5641 
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Chapter  5 

Verification  of  Technique 


5 . 1  300  mb  Independent  Data  Set 

The  basic  statistical  properties  of  the  variables 
included  in  the  123  cases  comprising  the  300  mb  independent 
data  set  are  shown  in  Table  10  for  comparison  with  the 
statistics  computed  for  the  300  mb  dependent  data  set  in 
Table  4. 

The  300  mb  SD  regression  equation,  including  the 
satellite  derived  parameters,  was  used  to  compute  estimates 
of  the  magnitudes  of  the  SD  extrema  of  the  independent  data 
set.  The  RMS  error  of  the  observed  versus  predicted  SD 
values  for  the  independent  data  set  was  33.27  meters.  A 
plot  of  the  residuals  for  the  independent  data  set  is  given 
in  Figure  13.  The  residuals  for  the  independent  data  set 
appear  to  have  no  trends.  Therefore,  the  300  mb  regression 
model  appears  to  be  adequate  and  stable  with  no  apparent 
systematic  error. 

5 . 2  1000/300  mb  Independent  Data  Set 

The  basic  statistical  properties  of  the  variables 
included  in  the  115  cases  comprising  the  1000/300  mb 
independent  data  set  are  shown  in  Table  11  for  comparison 
with  the  statistics  computed  for  the  1000/300  mb  dependent 
data  set  in  Table  7. 
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VAR 1ABIE 

MEAN 

STANDARD 

DEVIATION 

SKEWNESS 

KUR  TOS  I  S 

PLAT 

.  85865 

.  1  4973 

-.  18516 

- . 47763 

PLON 

. 70092 

. 36805 

.21312 

- . 79260 

SLAT 

.87734 

.  15327 

-  .  09850 

-  .  4627  1 

SLON 

.61980 

.36222 

.  19806 

- .58390 

AMP 

3 . 50483 

1 . 39591 

. 74888 

- . 24675 

RL 

2 . 88451 

1 .03516 

.77824 

- . 0764  1 

1ST 

2.34146 

.47420 

.66864 

-1 .55293 

SO 

•220.39840 

55.34015 

-  .08370 

.27860 

SR  U 

18.21138 

7.44219 

. 60404 

. 30371 

SR  V 

-  .  65854 

6.13378 

-.41914 

.246  1  2 

SVL 

3.55631 

3.91525 

-  .  70241 

-.21135 

SLH 

-75.08942 

91 .21101 

-.31323 

1 . 30044 

SLU 

5.52846 

6. 19116 

.01502 

. 35402 

SLV 

-.53659 

5.78095 

-  .  35753 

. 39805 

SLL 

24.97237 

20.31718 

.23392 

1.10163 

SVH 

71 . 13820 

179.50180 

-  .73296 

- .06124 

Svu 

12.89431 

3.3781 1 

.  19871 

- . 38882 

svv 

.21138 

.  65335 

1  .  32765 

4. 19123 

RLS 

2 . 1 7079 

.76418 

. 7607  ' 

. 32208 

R  L  5 

l  .67234 

.  29627 

.44464 

- .47865 

R  L  MA 

11.39112 

8 . 56090 

1  .  33408 

1.10512 

AMOR 

1 .21830 

.  23790 

. 77885 

1 .54160 

RLOA 

.85162 

. 16507 

.63477 

.654  1  5 

R  L  DC 

4 . 65501 

l .64013 

.62742 

. 15502 

RLS2 

1 . 66092 

.63788 

.66818 

. 42900 

123  TOTAL  OBSERVATIONS 
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PREDICTED  SD  VWJJE  (METERS) 


Figure  13.  Scatter  diagram  of  the  predicted  SD 
value  versus  the  residual  for  the 
300  mb  independent  data  set. 


VAR IABLE 


MEAN 


STANDARD 

DEVIATION 


SKEuNESS 


K  UR  TCS  I  S 


FLAT 

. 79690 

.15323 

-  .  20789 

.  3901  4 

FLON 

. 71623 

.  35677 

.23777 

-  .  97346 

SLAT 

.93593 

.  16644 

-  .  23297 

-  .  07636 

SLON 

.  59932 

. 34206 

.  18889 

-.99180 

AMP 

3.32952 

1.31369 

. 35556 

. 85699 

RL 

2 . 801 70 

.  92536 

.  48727 

-.61074 

1ST 

2 . 36522 

.48149 

. 55985 

-  1  .  68658 

SO 

-163.73910 

51 . 28595 

-  .04125 

-.51925 

SR  U 

16 . 10434 

6.25736 

. 55837 

.12241 

SR  V 

r .  26957 

6.11103 

- . 44373 

.  40937 

SVL 

2.44304 

3.93296 

-  .52790 

.  18706 

SLH 

-34.73912 

63.S38S4 

-.20703 

.  20020 

SLU 

3  .  18261 

4.91419 

.36437 

.03167 

SL  V 

-.24348 

5.51140 

-.38617 

.  26461 

SLL 

16 . 34540 

15.26734 

-.09495 

-  .05878 

S  VH 

103 . 30430 

157.83590 

- . 44509 

- . 36990 

svu 

1 3 . 20000 

3.65037 

. 16673 

- . 53047 

S  vv 

.33913 

.  74468 

1 .00293 

1.83109 

RLS 

2.04516 

.72214 

. 34224 

-  .  74936 

R  L  5 

1.65122 

.27413 

.21961 

- . 85480 

RLMA 

10.33154 

7.33653 

1.51171 

2.6672  1 

AMOR 

1  .  19500 

.28186 

1  .28466 

2 .426  10 

R  L  DA 

.8783  1 

.  1  8548 

. 1 7222 

- . 04934 

RLOC 

4.34616 

1.49327 

.31197 

- . 83532 

RLS2 

1  .52489 

.62007 

.  20735 

-  .  698  10 

115  TOTAL  OBSERVATIONS 


The  1000/300  mb  SD  regression  equation,  including  the 
satellite  derived  parameters,  was  used  to  compute  estimates 
of  the  magnitudes  of  the  SD  extrema  of  the  independent  data 
set.  The  RMS  error  of  the  observed  versus  predicted  was 
32.7  meters.  A  plot  of  the  residuals  for  the  independent 
data  set  is  given  in  Figure  14.  The  residuals  for  the 
independent  data  set  appear  to  have  no  trends.  Therefore, 
the  1000/300  mb  regression  model  appears  to  be  adequate  and 
stable  with  no  apparent  systematic  error. 
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Chapter  6 

Conclusions  and  Recommendations 

6 . 1  Conclusions 

The  primary  objective  of  the  present  study  was  to 
formulate  a  quasi-objective  technique  to  derive  the 
magnitudes  of  the  negative  extrema  of  the  300  mb  and 
1000/300  mb  height  and  thickness  short  wavelength  component 
fields  based  upon  satellite  observed  spiral  cloud  pattern 
geometry  and  various  long  wavelength  field  parameters.  It 
was  shown  that  quantitative  use  of  satellite  observed  spiral 
cloud  patterns  and  parameters  derived  from  the  long 
wavelength  fields  can  be  made  in  modifying  the  300  mb  and 
1000/300  mb  height  and  thickness  fields.  Proper 
interpretation  of  the  satellite  imagery  permits  the 
identification  and  location  of  the  short  wavelength  systems 
in  the  300  mb  and  1000/300  mb  fields.  The  regression 
equations,  including  the  satellite  parameters  as  predictors, 
permit  quantitative  estimates  of  the  magnitudes  of  the 
negative  extrema  in  the  short  wavelength  300  mb  and  1000/300 
mb  fields. 

The  multiple  regression  correlation  coefficient  for  the 
300  mb  data  shows  little  appreciable  difference  compared  to 
Nagle  and  Hayden's  correlation  coefficient  for  the  500  mb 
field  (0.31  versus  0.80  respectively).  The  root  mean 
squared  (RMR,  error  of  30.9  meters  for  the  500  mb 
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independent  data  set  is  appreciably  lower  than  the  RMS  error 
of  41  meters  for  Nagle  and  Hayden's  500  mb  independent  data 
set.  The  500  mb  height  analysis  using  surface,  radiosonde 
height  and  satellite  vertical  temperature  profiles  has  a  RMS 
error  of  10.6  meters  when  computed  for  the  relatively  dense 
data  network  of  Northern  Europe  (Murphy  and  Williamson, 

1976).  This  can  be  considered  as  the  inherent  error  in 
the  current  "perfect"  analysis.  The  RMS  error  value  of  the 
predicted  (the  "first  guess")  300  mb  heights  as  given  by 
Rutherford  (1976)  is  32.7  meters  computed  for  the  entire 
northern  hemisphere.  This  implies  that  the  RMS  height  error 
over  data  sparse  areas  can  exceed  40  meters.  Therefore,  the 
technique  developed  in  the  present  study  can  significantly 
improve  the  accuracy  of  the  300  mb  geopotential  height  field 
in  data  sparse  regions. 

Upon  examining  the  1000/300  mb  thickness  field  RMS 
height  error  of  32.7  meters  in  light  of  the  above 
discussion,  it  appears  that  the  present  technique  can  also 
significantly  aid  the  analysis  in  data  sparse  areas. 

Equally  important  is  that  the  1000/300  mb  thicknesses 
estimated  by  this  technique  can  insure  vertical  consistency 
of  the  mass  field  which  is  an  essential  requirement  for 
initializing  a  multi-level  primative  equation  model. 

The  regression  equations  developed  excluding  the 
satellite  parameters  can  be  used  when  specification  of  the 
cloud  parameters  is  ambiguous.  The  estimates  of  the  3C 


values  computed  with  these  regression  models  are 
conservative  estimates  with  larger  RMS  errors  than 
associated  with  the  regression  models  including  the 
satellite  parameters.  These  equations  could  be  used  for 
systems  in  the  formative  stage  of  development  when  the 
interpretation  of  the  satellite  image  is  difficult,  and 
requires  considerable  subjectivity. 

Therefore,  provided  the  satellite  images  are  properly 
interpreted  and  that  the  long  wavelength  component  fields 
provided  by  the  12-hr  numerical  forecasts  can  be  accepted  as 
a  reasonable  base,  the  technique  developed  in  the  present 
study  can  significantly  improve  the  accuracy  of  the  300  mb 
and  1000/300  mb  geopotential  height  and  thickness  fields  in 
data  sparse  regions.  Appendix  C  provides  the  details  for 
adapting  these  techniques  for  operational  usage. 

6 . 2  Recommendations  for  Future  Study 

Verification  of  the  zero  line  positioning  and  testing 
of  the  techniques  to  define  the  spatial  characteristics  of 
the  300  mb  and  1000/300  mb  SD  fields  should  be  further 
explored.  This  will  require  satellite  imagery  [preferably 
from  a  geostationary  satellite,  in  order  to  avoid  problems 
associated  with  mosaics)  valid  at  analysis  times.  The 
spatial  errors  between  the  satellite  image  and  the  field 
analysis  would  then  be  at  a  minimum. 

The  impact  of  this  technique  on  a  regional  or  global 
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numerical  model,  and  the  ability  of  the  objective  analysis 
and  initialization  schemes  to  retain  the  adjustments  to  the 
fields  should  be  pursued.  However,  if  the  impact  is  neutral 
or  negative,  an  explanation  may  be  that  the  additional 
satellite  information  provided  by  this  technique  are  not 
being  exploited  fully  by  the  existing  system.  This  is  quite 
probable  if  the  existing  system  has  been  designed  primarily 
for  assimilating  synoptic  rawinsonde  data.  Finally, 
assessments  of  the  impact  should  be  qualitative  as  well  as 
statistical,  should  discriminate  according  to  geographical 
area  and  space  time  scales,  and  should  give  emphasis  to 
synoptically  active  test  situations. 


APPENDIX  A 


where  <p  is  the  smoothed  field,  is  the  original 

2 

field,  <5  is  a  constant,  7  the  two-dimensional 
Laplacian  operator,  and  a  is  the  degree  of  smoothing 
( Holl ,  1963),  Equation  (Al)  may  be  expressed  in  a  more 
general  operator  notation  as 


0  [a]  =  <p  [  0  ]  -  f>  [  0  ,  ot  ] 


( A2 ) 


where  f[a]  is  the  field  ;  smoothed  to  degree  of 
smoothing  a  ,  and  $[0,a]  is  the  residual  field 
resulting  from  the  smoothing  operating  on  the  :  field 
from  its  initial  state  a  =  0  to  the  final  state 
a  =  a  . 

The  properties  of  the  smoother  can  be  examined  by 
assuming  the  $  field  can  be  represented  by  a  one¬ 
dimensional  Fourier  component  with  wave-number  k  and 


whose  amplitude  Ak  will  be  a  function  of  the  degree 
of  smoothing.  Equation  (Al)  can  be  written  in  differen¬ 
tial  form  as 


£±  =  6  2  7  2 

3a 


The  one-dimensional  Fourier  component  is  expressed  by 


MO]  =  Ak[0]  e' 


where  i  =  /-T  .  Substituting  (A4)  into  (A3)  and 
solving  yields 


=  MO]  e 


.2-2 
-k  3  a 


Ak [ a ]  =  Ak [ 0 ]  e 


-k2  5  2a 


The  amplitude  reduction  factor  which  is  a  function 
of  the  degree  of  smoothing  and  wave-number  is  given  by 


R  = 


Ak[M 

VO] 


.2  .2 
-k  i  a 


The  larger  the  degree  of  smoothing  (large  a)  the  more 
the  amplitude  will  be  reduced.  It  is  readily  seen  that 
for  a  given  degree  of  smoothing  the  shorter  waves' 


(large  k)  amplitudes  are  greatly  reduced. 


2.  Analysis  of  the  Numerical  Analogue 


Computationally  Equation  (Al)  is  integrated  by 

repeated  applications  of  the  Laplacian  operator  with 

varying  increments  in  a  (  a.  ,  a.  ,  .  .  .  , a  )  .  Using 

l  Z  n 

Equation  (A5)  one  can  show  that  the  resultant  smoothing 
is  cumulative  in  a  /  since 


v2  -2 

$[aL]  [a2]  =  $[0]  [a2]  e“*  0  31 


=  „[0]  e-^2 ( al  +  a2 } 


=  $ +  a2 ] 


The  explicit  finite  difference  analogue  of  Equation  (A8) 


=  2(0]  +  *  : [ 0 ] 


( +  a2  ] 


^  V  (  i  ^  ^  2  * 


or  in  general  for  the  nth  degree  of  smoothing: 


'  ;  IZ_  J1  hn 


where  is  the  five-point  finite  difference  approxi¬ 

mation  of  the  Laplacian  operator  on.  a  lattice  with 
gndlength  5  .  For  reasons  of  numerical  stability  the 
a  Vs  increase  with  successive  steps. 

The  smoothing  of  degree  a  is  complete  when  the 
sum  of  the  a  Vs  reaches  the  predetermined  value, 
cp  [ ot ]  is  the  long-wave  component  field  and  p  [ 0 j  -  y  [a] 
is  the  short-wave  component  field. 

Equation  (A7)  can  be  used  to  construct  a  filter 
(smoother)  response  curve  where  the  wave-number  is 
replaced  by  wavelengths  in  grid  increments,  i.e. 

k  =  2^(n6)_1  .  (All) 


The  amplitude  reduction  factor  becomes 


R 


a/n 


2 


(A12) 


Using  Equation  (A12)  a  family  of  curves  corresponding 
to  different  a  for  the  reduction  factor  as  a  function 
of  n  can  be  constructed.  However,  a  single  more 
general  curve  can  be  derived  by  introducing  a  new 


variable  x  which  is  a  function  of  both 


X 


and  n  : 


where  r  is  an  arbitrary  positive  constant  used  to 
scale  x  .  Equation  (A12)  can  be  rewritten  as 


For  plotting  purposes  the  constant  r  was  chosen  by 
requiring  that  R  =  0.5  when  x  =  1.0  ,  which  leads 
to  the  numerical  value  for  r  of  0.693147  .  The 
filter  curve  for  this  scheme  is  shown  by  the  solid  line 
in  Figure  15  . 

An  ideal  filter  would  be  a  vertical  line  at  x  =  1.0 
The  filtering  technique  described  above  is  known  as 
first-order  smoothing  and  is  not  as  sharp  as  the  ideal 
filter.  The  resulting  long-wavelength  field  contains 
some  short-wavelength  components  and  vice  versa.  The 
application  of  higher  order  filters  can  compensate  this 
undesirable  feature  at  the  expense  of  computation  time. 

The  concept  of  higher  order  smoothing  is  demonstrated 
with  the  following  example.  Suppose  that  after  applying 
first-order  smoothing,  the  short-wavelength  field 
contains  95  percent  of  the  amplitude  of  a  given  short 
wave  and  erroneously  contains  15  percent  of  the  amplitude 
of  a  given  long  wave,  due  to  the  imperfect  first  order 


RRF 


EXP ( -R*X*X) 


If  the  filtering  process  is  now  applied  to  the  short- 
wavelength  field,  to  yield  large  and  small  scale  components, 
the  amplitude  of  the  long  wave  will  be  reduced  in  the  second 
short-wavelength  component  field  to  0.15  x  0.15  or  2.25 
percent  of  its  original  amplitude.  The  short  wave  amplitude 
will  be  reduced  in  the  second  short-wavelength  component 
field  to  0.95  x  0.95  or  90.25  percent  of  its  original 
amplitude.  This  second-order  smoothing  is  much  sharper  than 
the  first  order  because  the  longer  wavelengths  are  filtered 
more  rapidly  than  the  shorter  wavelengths. 

3.  Second  Order  Smoothing 

Equation  (A2)  for  first-order  smoothing  can  be 
rewritten  as: 

y[0»a]  =  <P(0]  -  MM  •  ( A 1 5 ) 

Second-order  smoothing  involves  smoothing  the  residual 
field  a  second  time: 

5  [  0  ,  a  ]  [0,a]  =  ‘M0]  -  Mi]  ;  [0,  l] 

=  MO]  [0]  -  M0]  [a]  -  MM  [0]  +  M  M  ( M. 

(A16) 


60 


From  Equation  (AS) 


;[0,ij(0,i]  *  o  ( 0  ]  -  2  $  [  a  ]  +  ; [ 2a ] 

Adopting  the  convention 

0  (a)  [aj  *  $  ( 2a)  =  o  [a]  2. 

Equation  (A17)  can  be  written  as: 


( A17 


( A18 


*C0,a]2  ,  *(i  -  [a] )  2  (j 

or  in  general 

♦  [0,a]n  «  o(i  -  £a])n  >  (p 

The  resulting  long-wavelength  component  field  after  nth 
degree  of  smoothing  is: 


♦  LW  ■  ♦CO]  -  ♦ i 0  ,  a  ]  n 

-  ill  -  (1  -  C*])n)  .  (A 

The  smoothing  process  used  for  this  study  was  third- 
order  smoothing.  The  long-wavelength  field  is  given  by 

5  w  »  3ila)  -  3$ I  2a)  +  i 1 3a)  .  (A22 


Equation  (A22)  shows  that  in  practice  the  field  is 
smoothed  to  3  a  rather  than  by  repeated  smoothing, 
which  computationally  xs  more  efficient. 

The  amplitude  reduction  factor  for  third-order 
smoothing  is 


-  3e 


-  3e 


-  2rx 


-  3rx‘ 


By  requiring  that  R  =  0.5  when  x  =  1.0  the  numerical 
value  of  the  constant  r  is  1.5784264  .  Figure  15 
shows  that  third  order  smoothing  is  indeed  sharper  than 
first  order.  Using  Equation  (A23)  and  the  definition 
of  r  ,  the  corresponding  degree  of  smoothing  is 
a  =  8.52  . 

Equation  (Al)  defines  the  procedure  for  explicit 
smoothing.  However,  implicit  smoothing  is  found  to  be 
more  efficient  and  is  used  m  practice.  One  smoothing 
cycle  is  represented  by 


Equation  (A24)  is  solved  by  standard  over relaxa  t  icr. . 
The  optimum  relaxation  factor  (Holl,  1967)  *  .3 

given  by : 


where 


q  =  2  5  a/ (25a  + 1)  .  (A26 

The  number  of  iterations  required  for  adequate  accuracy 
is 


N 


in  0.02 

-In  (ui  -  1) 


(A27 


The  6a  for  each  smoothing  cycle  is  defined  such  that 

after  15  cycles  the  fields  will  be  smoothed  to  a  ,  2a, 

and  3a  .  A  normalized  table  for  5a  is  given  in 

* 

Table  12  .  The  actual  5a  used  for  each  cycle  is 

found  by  multiplication  of  the  normalized  value  and 
the  desired  degree  of  smoothing. 


TABLE  12 


Normalized  smoothing  increments  (after 


Cycle  No. 


0 . 006 

: . '  "6 

0 . 009 

■'  •  ~  L  5 

0.014 

"  .  ^29 

0.020 

0.049 

0.030 

0 . 079 

0.045 

0.124 

0 . 068 

0 . 192 

0.102 

0 . 294 

0.154 

0.448 

0.230 

0.678 

0.322 

1.000 

0.500 

1  .  500 

0.500 

2 . 000 

0 . 500 

2 .500 

0.500 

3 . 000 

APPENDIX  B 


The  linear  regression  of  y  on  the  variable  x  can  be 
extended  to  the  multiple  regression 


=  BiXl  +  62X2  +  ...  +  -nXn  =  S 


*3X3  - 


(  B 1  ) 


3  =  1 


where  X),X^,..*,Xn  are  the  n  different  variables  and  £j, 
02  /•••»  6n  are  the  model  parameters  of  the  regression  line 
The  special  case  in  which  Xi  is  identically  one,  would 
mean  that  8j  had  now  assumed  the  role  of  the  regression 
constant  or  y-intercept.  Suppose  that  M  experimental 


observations  are  made  and  that  X  j  j  is  the  value  of  the 


variable  X-  at  the  ith  observation.  Letting  Yj  be  the 


ith  observed  value  of  y  ,  the  formulation  of  the  multiple 
regression  procedure  requires  one  to  minimize  the  sum  of 
squares 

M  N 

S  =  J  <*i  -  E  a:Xij)2  •  (32) 


i=l  3=1 

Setting  3s/^Sk  =  0  for  '<  =  l,2,...,n  ,  and  replacing 
the  S j  ,  with  their  least-squares  estimates  b  ,  we  have 


M 


n 


I 


X i k  (Yi  -  bjXij)  =  0 

i=l  3=1 


That  is, 

n  M  M 

£bi  Z  xikYiJ  ■  £ 

:=i  1=1  1= l 

k  =  1,2 . n  . 

Now  define  X  as  the  Mx.n  matrix 

its  ith  row  and  jth  column: 


X  —  (  x  i  j  ) 


r  k  foil  jW  u:q  column  vectors  are  also 

y  =  CY1 

b  =  ’  d  ^ 

V  he:  ■  -  he  superscript  t  ienmes  the 
/ector.  The  system  of  normal  eauations 


Y2 
n  -> 


X  i  <  Y  j 


IK1!  < 


con 


taming  x 


XM,  n 


define  1 : 


Yn 


1  - 


n  n 


i  n  s  no  s  * 


(  B  3  )  :  ~i  n  h 


which  has  the  solution 


b  =  (Xt  x)"1  xt  y  . 

The  Mxn  symmetrical  var iance-covar iance  matrix  V  is 
defined  by 


V  »  (vij)  =  (*t  X)-l  a2  ,  ( B9 ) 

where  ^2  is  estimated  from 

Yt  y  -  bt  Xt  y 

s2  a:  M  -  n  •  (BIO) 

The  diagonal  and  off-diagonal  elements  of  V  give, 

respectively,  the  variances  and  covariances  of  the  bj  * 
that  is, 

var(bj)  =  v i j  ,  and  cov(bi,bj)  =  vij  . 

The  system  of  normal  equations  ( B 3 )  can  be  solved  by  any 
number  of  conventional  methods  such  as  Gaussian  elimination 
or  matrix  inversion. 

In  the  present  study,  the  Biomedical  Computer  Programs 
P-Series  (BMDP-79)  library  developed  at  the  Health  Sciences 
Computing  Facility  at  the  University  of  California  at 


67 


Los  Angeles  was  used  (Dixon  and  Brown,  1979).  The  stepwise 
multiple  regression  program  P2R  was  employed  for  the 
regression  analysis.. 

P2R  computes  estimates  of  the  parameters  of  a  multiple 
linear  regression  equation  in  a  stepwise  manner.  The 
variables  are  entered  (forward  stepDing)  or  removed 
(backward  stepping)  from  the  equation  one  at  a  time 
according  to  any  of  four  possible  criteria.  Univariate 
statistics  computed  for  each  variable  are:  mean,  standard 
deviation,  coefficient  of  variation,  skewness,  kurtosis, 
maximum,  and  minimum  values.  Results  are  printed  at  each 
step  along  with  the  multiple  correlation  coefficient  R 
(i.e.  the  correlation  of  the  dependent  variable  y  with  th< 
predicted  value  y  )  as  well  as  the  multiple  R2  ,  the 
adjusted  Ra2  =liven  bY 


,  * 


rA2  =  R2  -  n ( 1  -  R2 ) / ( M  -  n)  (311) 

(Thiel,  1971),  and  the  standard  error  estimate  given  by: 


M 


The  analysis  of  variance  table  containing  the  regression 
sums  of  squares,  residual  sums  of  squares  an  1  the  F 


(>  8 


I*.  •' .  •  ,  •’  .  *  ,  *  .  *  » 

—  .m  -V  ~  l.U» 
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inclusion  in  the  equation  in  a  manner  that  maximizes  the  P 
ratio,  the  level  of  significance  of  the  P  ratio  cannot  be 
obtained  from  an  P  distribution. 

Any  one  of  four  methods  can  be  used  for  entering  or 
removing  variables  from  the  equation  at  each  step.  The 
methods  are: 


P  The  variable  with  the  smallest  F-to-Remove, 

given  by: 

SS  (residuals  if  variable  is  removed)  -  SS  (residuals) 

SS  (residuals)  )  (M  -  n) 

Where  SS  is  the  sum  of  squares,  is  removed  it 
its  F-to-remove  is  lower  than  the  specified  F- 
to-remove  limit.  If  no  variable  meets  this 
criterion,  the  variable  with  the  largest  F-to- 
enter  defined  by 

SS  (residuals)  -  SS  (residuals  after  next  step)  , 

SS  (residuals  after  next  step) /(M-n-1 ) 

is  entered  if  the  P-to-enter  value  exceeds  the 
specified  P-to-enter  limit. 

SWAP  Same  as  method  P  above,  except  when  no  variable 
meets  the  F-to-remove  criterion,  a  variable  in 
the  equation  is  exchanged  with  a  variable  not 


the  multiple  regression  correlation  coefficient 
R.  If  no  variable  can  be  exchanged,  a  variable 
is  entered  as  in  P. 

R  The  variable  with  the  smallest  F-to-remove  is 

removed  if  its  removal  results  in  a  larger 
multiple  R  than  was  previously  obtained  for 
the  same  number  of  variables.  If  no  variable 
meets  this  criterion,  a  variable  is  entered  as 
in  F. 

RSWAP  Same  as  method  R,  except  when  no  variable  meets 
this  criterion,  a  variable  in  the  equation  is 
exchanged  with  a  variable  not  yet  in  the 
equation  if  the  exchange  increases  the  multiple 
correlation  coefficient  R.  If  no  variable  can 
be  exchanged,  a  variable  is  entered  as  in  F. 

The  method  F  requires  the  fewest  computations. 

Since  all  four  methods  begin  with  no  variables  in  the 
equation,  they  appear  to  describe  a  forward  stepping 
algorithm.  The  F-to-enter  and  F-to-remove  limits  "'an  e a ■> h 
be  assigned  two  values.  The  first  pair  is  used  until  all 
variables  that  meet  the  criteria  have  been  entered;  th° 
second  pair  of  limits  is  then  use  i  to  rem  v:o  tne  variable 
(backward  stepping'.  I?y  assigning  very  low  values  to  ‘ h 


PROGRAM  CONTROL  INFORMATION 


/PROBLEM 

/INPUT 

/VARIABLE 

/TRAN 


/REGRESS 


/PRINT 

/END. 


TITLE  IS  300  MB  SD  SINAP  REGRESSION  . 

VARIABLES  =  19. 

FORMAT  =  (  6X,6E14 . 6, 2X,  II  ,  /  ,  i  2  F.  10. 3  )  . 

UNIT  =  20. 

NAMES  ARE  F L AT , F LON, SL AT, SLON, AMP, RL , I S T , SD,  SRU  , 
SRV,SVL,SLH,SLU,SLV,SLL,SVH,SVU,SVV,SVL,RLS,RL5,  j 
RLMA, AMDR,RLDA,RLDC,RLS2.  1 

ADD  =  7. 

RLS  =  RL*S I N ( SLAT )  . 

RL  5  =  SORT(PL). 

RLMA  =  R  L  *AMP  . 

AMOR  =  AMP/RL. 

R  L  DA  =  RL/AMP. 

RLDC  =  RL/COS( SLAT) . 

RLS2  =  RL*S  I  N  ( SLAT ) *S  I  N ( SLAT) . 

DEPEND  =  SD. 

INDEP  =  SLAT,  SLON,  AMP,  RL, SRU, SRV, SVL, SLH, SLU, SLV, 
SLL,SVH,SVU,SVV,SVL, RLS, RL5,FLMA,AMDR,RLDA, RLDC, 
RLS2. 

TITLE  -  FORWARD  STEPPING  MULTIPLE  REGRESSION  . 
METHOD  =  F. 

ENTER  =  4.0. 

REMOVE  =  3.9. 

TOL  =  0.05. 

DATA. 

CORR. 


BM DP  UNIT  NO.  20  SPECIFIED  IN  THE  INPUT  PARAGRAPH  WILL  REFER  TO 
LOCAL  FILE  NAME  TAPE20  FOR  THIS  PROBLEM. 


Fiuure  16.  P2R  oroaram  control  Lmmiaa". 


APPENDIX  C 


The  operational  application  of  the  technique  developed 
in  the  present  study  and  described  herein,  is  designed  to 
utilize  the  Naval  Environmental  Prediction  Research 
Facility's  Satellite-Data  Processing  and  Display  System's 
(SPADS)  unique  capabilities  of  allowing  interactive  digital 
image  processing  (Schramm,  Zeleny,  Nagle,  and  Weinstein, 
1982).  Graphic  overlays  of  the  meteorological  fields 
superimposed  upon  the  satellite  image  make  this  system  ideal 
for  the  application  of  the  present  technique. 

Since  parameters  derived  from  the  long  wavelength  field 
are  required  for  the  regression  models  for  estimating  the 
magnitudes  of  the  negative  extrema  of  the  SD  field,  it  must 
be  assumed  that  the  long  wavelength  field  is  accurately 
specified.  In  a  realtime  application  the  current  long 
wavelength  analysis  is  not  availiable,  but  rather  the  12-hr 
numerical  forecast  of  the  fields  are  to  be  used.  The 
assumption  being  that  the  12-hr  forecast  .(the  "first  guess") 
contains  sufficient  information  for  accurate  specification 
of  the  long  wavelength  fields.  Nagle  and  Hayden  tested  this 
assumption  by  comparing  the  differences  between  the  12-hr 
first  guess  field  and  subsequent  analysis  for  500  mb  height 
fields.  They  also  compared  tne  differences  of  the 
Laplacians  of  the  fields  and  found  in  both  cases  the  errors 
were  acceptably  low  for  use  in  the  current  context. 


The  primary  satellite  data  used  will  be  digital 
satellite  imagery  obtained  from  the  geostationary  satellite 


of  the  GOES  series.  The  12-hr  numerical  300  mb  and  1 000/31 
mb  forecast  fields  are  computer  processed  to  produce  the  St 
SR,  SL  and  SD  fields.  Selection  of  the  GOES  image  sector 
within  the  hour  preceeding  the  analysis  time  (i.e.  1200  or 
000  GMT)  is  then  made  and  the  image  displayed.  The  analyst 
then  interprets  the  spiral  cloud  system  by  defining  the 
parameters  required  by  the  regression  equation.  Once  the 
analyst  is  satisfied  with  the  defined  spiral  parameters,  th 
value  of  the  3D  extrema  is  then  computed. 

The  zero  line  defining  the  region  of  influence  of  the 
computed  3D  feature  is  then  defined  along  the  amplitude  and 
wavelength  axes  designated  by  graphical  means  on  the  color 
monitor.  The  zero  line  intersects  these  lines  at  the  inner 
edge  of  the  major  frontal  cloud  band.  This  defines  the 
major  and  minor  axes  of  an  ellipse  that  can  b<=  used  to 
approximate  the  3D  feature.  Another  possibility,  would  b^ 
to  simply  define  the  position  of  the  zero  line  with  the 
graph  pen  and  pad.  The  3D  feature  coul  i  "her.  be  define!  by 
“he  following  relationship 


-  L 


mm 
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where  Z  is  the  value  at  the  grid  point  defined  by  r,  where 
is  the  distance  from  the  center  to  “hQ  gr :  d  p-in",  ?.  . 
distance  from  the  center  to  the  zero  line  passing  “hr : ugh 
“he  grid  point  and  is  “  rm  ."j  v-alu-  --~.ru-  :  ■; 
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regression  equation.  Figure  17  is  a  schematic 
representation  of  the  above. 

Once  the  region  of  influence  is  defined,  the  SD  values 
at  the  grid  points  within  the  region  of  influence  when  added 
to  the  SR  field  compose  an  array  of  bogus  height  or 
thickness  values  to  be  blended  with  the  rawinsonde  and 
aircraft  data  used  for  the  analysis.  This  procedure  can 
also  be  used  for  the  6-hr  upd.ate  cycles  when  an  explosive 
case  of  cyclogenesis  is  apparent  in  the  satellite  imagery 
but  not  reflected  in  the  twice  daily  analysis  times.  A 
schematic  of  this  SINAP^-  process  is  shown  in  figure  18. 


1  SINAP: 
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Figure  18. 


Schematic  flow  diagram  for  the 
operational  usage  of  the  SINAP 
techn ique . 
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SILVER  SPRING,  MD  20910 


COMMANDER 
PACMISTESTCEN 
GEOPHYSICS  OFFICER 
PT.  MUGU,  CA  93042 


5WW/DN 

LANGLEY  AFB,  VA  23665 


COMMANDER  &  DIRECTOR 
ATTN:  DELAS-D 
U.S.  ARMY  ATMOS.  SCI.  LAB 
WHITE  SANDS  MISSILE  RANGE 
WHITE  SANDS,  NM  88002 


DIRECTOR 

NATIONAL  METEORO.  CENTER 
NWS,  NOAA 
WWB  W32,  RM  204 
WASHINGTON,  OC  20233 


DIRECTOR 

PACIFIC  MARINE  CENTER 
NATIONAL  OCEAN  SURVEY,  NOAA 
1801  FAIRVIEW  AVE.,  EAST 
SEATTLE,  WA  98102 


CHIEF 

MESOSCALE  APPLICATIONS  BRANCH 
NATIONAL  EARTH  SAT.  SERV. 

1225  W.  DAYTON 
MADISON,  WI  53562 


HEAO,  ATMOS.  SCIENCES  DIV. 
NATIONAL  SCIENCE  FOUNDATION 
1800  G  STREET,  NW 
WASHINGTON,  OC  20550 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
ATTN:  CODE  44 
OAHLGREN,  VA  22448 


USAFETAC/TS 
SCOTT  AFB,  IL  62225 


OFFICE  OF  STAFF  METEOROLOGY 
WESTERN  SPACE  &  MISSILE 
CENTER  (WE) 

VANDENBERG  AFB,  CA  93437 


DIRECTOR  (12) 

DEFENSE  TECH.  INFORMATION 
CENTER,  CAMERON  STATION 
ALEXANDRIA,  VA  22314 


DIRECTOR 

NATIONAL  EARTH  SAT.  SERV/SEL 
FB-4,  S321B 
SUITLAND,  MD  20233 


NOAA-NESDIS  LIAISON,  CODE  SC2 
NASA-JOHNSON  SPACE  CENTER 
HOUSTON,  TX 


MANAGER 

NATIONAL  EARTH  SAT. SERV. 
SATELLITE  FIELD  SERV.  STN. 
701  C  ST.,  BOX  45 
ANCHORAGE,  AK  99513 


COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN:  DR.  WILLIAM  GRAY 
FORT  COLLINS,  CO  80523 


NAVAL  SPACE  SYSTEMS  ACTIVITY 
CODE  60 

P.O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 


AFGL/LY 

HANSCOM  AFB,  MA  01731 


OFFICER  IN  CHARGE 
SERVICE  SCHOOL  COMMAND 
DET.  CHANUTE/STOP  62 
CHANUTE  AFB,  IL  61868 


DIRECTOR,  SYSTEMS  DEVELOPMENT 

wuc  NDAA 

RM  1216  -  THE  GRAMAX  BLDG. 
8060  13TH  ST. 

SILVER  SPRING,  MD  20910 


DIRECTOR 

NATIONAL  SEVERE  STORMS  LAB 
1313  HALLEY  CIRCLE 
NORMAN,  OK  73069 


METEOROLOGIST  IN  CHARGE 
WEA.  SERV.  FCST.  OFFICE,  NOAA 
660  PRICE  AVE. 

REDWOOD  CITY,  CA  94063 


DIRECTOR 

TECHNIQUES  DEVELOPMENT  LAB 
GRAMAX  BLDG. 

8060  1 3TH  ST. 

SILVER  SPRING,  MD  20910 


CHAIRMAN,  METEOROLOGY  DEPT. 
CALIFORNIA  STATE  UNIVERSITY 
SAN  JOSE,  CA  95192 
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COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN:  LIBRARIAN 
FT.  COLLINS,  CO  80523 


DIRECTOR,  REMOTE  SENSING  LAB 
UNIVERSITY  OF  MIAMI 
P.O.  BOX  248003 
CORAL  GABLES,  FL  33124 


CHAIRMAN,  METEOROLOGY  DEPT. 

UNIVERSITY  OF  UTAH 

SALT  LAKE  CITY,  UT  84112 


LIBRARIAN,  METEOROLOGY  DEPT. 
UNIVERSITY  OF  MELBOURNE 
PARKVILLE,  VICTORIA  3052 
AUSTRALIA 


CHAIRMAN,  METEOROLOGY  DEPT. 
METEORO.  &  SPACE  SCI.  BLDG. 
1225  W.  DAYTON  ST. 

MADISON,  HI  53706 


CHAIR, MAN,  METEOROLOGY  DEPT. 
PENNSYLVANIA  STATE  UNI V. 

503  DEIKE  BLDG. 

UNIVERSITY  PARK,  PA  16802 


SCIENCE  APPLICATIONS,  INC. 
2999  MONTEREY-SALINAS  HWY. 
MONTEREY,  CA  93940 


BUREAU  OF  METEOROLOGY 
BOX  1289K,  GPO 
MELBOURNE,  VIC,  3001 
AUSTRALIA 


UNIVERSITY  OF  WASHINGTON 
ATMOSPHERIC  SCIENCES  DEPT. 
SEATTLE,  WA  98195 


CHAIRMAN 

METEOROLOGY  DEPT. 
MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 
CAMBRIDGE,  MA  02139 


MR.  W.  G.  SCHRAMM/WWW 
WORLD  METEOROLOGICAL 
ORGANIZATION 

CASE  POSTALE  #5,  CH-1211 
GENEVA,  SWITZERLAND 


INSTITUTE  OF  OCEANOGRAPHY 
UNIV.  OF  BRITISH  COLUMBIA 
VANCOUVER  BC,  CANADA  V6T-1W5 


EUROPEAN  CENTRE  FOR  MEDIUM 
RANGE  WEATHER  FORECASTS 
SHINFIELD  PARK,  READING 
BERKSHIRE  RG29AX,  ENGLAND 


EUROPEAN  SPACE  OPERATIONS  HEAD,  DATA  PROCESSING  SEC. 

ATTN:  DR.  J.  MORGAN,  METEO.  GERMAN  MILITARY  GEOPHYS. 
SAT.  DATA,  MANAGEMENT  DEPT.  MONT-ROYAL,  D5580 
R.  BOSCH  STR  5  D61  DARMSTADT  TRAVEN-TRARBACH 
FEDERAL  REPUBLIC  OF  GERMANY  FEDERAL  REPUBLIC  OF  GERMANY 


